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We present preliminary results for the D s meson spectrum and decay constants in unquenched lattice QCD. 
Simulations are carried out with 2+1 dynamical quarks using gauge configurations generated by the MILC 
collaboration. We use the "asqtad" a 2 improved staggered action for the light quarks, and the clover heavy quark 
action with the Fermilab interpretation. We compare our spectrum results with the newly discovered + and 1 + 
states in the D s system. 



1. INTRODUCTION 

D physics is assuming a larger importance in 
lattice QCD and in CKM phenomenology with 
the arrival of the CLEO-c charm factory pQ. 
CLEO-c will measure the leptonic decay con- 
stants fn and fr>, to an accuracy of around 2%. 
It will measure the amplitudes of the semileptonic 
decays D — ► Ttlu and D — > Klv to an accuracy 
of around 1%. This will produce new determina- 
tions of the CKM matrix elements V c d and V cs to 
the accuracy that can be achieved in the required 
lattice calculations, and new checks of the uni- 
tarity triangle. New, precise tests of lattice QCD 
will come from the amplitude ratios Jd/D — > izlv 
and JdJD — > Klv. These ratios, which are in- 
dependent of the CKM matrix, will provide the 
most precise tests in existence of the types of lat- 
tice QCD calculations required to extract CKM 
matrix elements from B physics and D physics. 

The recent discovery of the positive parity part- 
ners of the D s and the D* [2] has added new in- 
terest to the spectrum of the D s system. The new 
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states lie significantly below quark model predic- 
tions, and below the DK threshold, so the states 
are quite narrow. 

2. METHODS 

We use the "asqtad" order a 2 improved stag- 
gered light quark action, and the order a im- 
proved Fermilab heavy quark action. In un- 
quenched calculations, chiral extrapolation is the 
least well controlled remaining source of error, 
staggered fermions are the only method currently 
able to reach the region mi ~ m s /5 that seems 
to be required to control this error. Using un- 
quenched calculations that reach this region, a 
number of simple quantities that disagree with 
each other at the 10% level in the quenched ap- 
proximation come into good agreement We 
use the public MILC unquenched configurations, 
with two light and one strange sea quarks, with 
properties 

• 20 3 x 64, a -1/8 fm, 

• (m;,m s ) = (0.007,0.05), (0.01,0.05), 
(0.02,0.05), (0.03,0.05), (True m s « 0.041.), 
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• ~ 500 configurations at each mass, 4 time 
sources each, 



as described in 0]. 



Figure 2. the D s correlation function with the 
leading exponential scaled out. The exponential 
decay of the oscillating signal gives the mass of 
the + state. 



Figure 1. The D s correlation function. 



With staggered and naive light quarks, the + 
and 1 + are in the D s and D* correlators automat- 
ically and unavoidably because of fermion dou- 
bling. Naive fermions possess a doubling sym- 
metry, -0 - > (*7m75)( — that implies that 
quarks come in multiples with identical prop- 
erties. Local fermion operators connect to all 
of these doubled modes. Therefore, because of 
the doubling symmetry, an pseudoscalar current 
couples to a positive parity scalar state as well: 

^hlsip -> */i75(«7o75)(-l)V = i^hloiPi-lY' 
Oscillating signals for the positive parity states 
are present in the correlation functions for the 
D s and D* mesons, but they are small and hard 
to see, as can be seen in fig. ^ Fig- 121 shows the 
same correlation function with the leading expo- 
nential scaled out. The decay of the oscillating 
signal gives the mass of the + state. Because 
it is an excited state, the presence of a plateau 
is more difficult to ascertain than for a ground 
state. 



3. SPECTRUM 

The D s spectrum on the mi,m s = 0.01,0.05 
lattices is shown in fig. [3] Very little dependence 
on the sea quark mass is observed. The + and 
1 + states lie somewhat above experiment. How- 
ever, they are excited states and are somewhat 
more sensitive to the choices of the priors for 
the higher states in the Bayes fits than are the 
ground states. We obtain 0.86 (2) for the hyper- 
fine splitting divided by its experimental value. 
This may be compared with 0.82 (2) for the hy- 
perfine splitting in charmonium. The most likely 
source for these errors is the one-loop correction 
to 0(a) ipTi ■ Bijj operator in the clover action. 
Since this operator contributes twice to the char- 
monium hyperfine splitting and once to the the 
D s hyperfine splitting, it is reasonable to expect 
a larger effect in charmonium. 

4. LEPTONIC DECAY AMPLITUDE 

This work is part of a larger project with the 
MILC collaboration. Preliminary results for the 
amplitude for leptonic decay of the D s are shown 
in fig.0]as a function of the sea quark mass (using 
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Figure 3. The D s spectrum on the m;,m s 
0.01,0.05 lattices. 
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m s — 0.041). There is little dependence on the 
sea quark mass, as expected. The one-loop Fer- 
milab heavy-staggered light axial current renor- 
malization is in progress To make an esti- 
mate of the one-loop correction, we use the for- 
mula Z% 1 = pas/z^zJ 0, using Z$ h = 1.33(2) 

and Z$ = 0.86(5) [7 . This leads to a current 
result for the decay constant of fr>, = 240 MeV 
+/- 0(a). However the order a correction could 
be large, potentially as large as 30%. This is to 
be compared with Ryan's unqucnched world av- 
erage of 250 (30) MeV [8] and the unquenchcd 
NRQCD heavy staggered light result of 289(a 2 ) 
0. The CKM independent quantity f D Jf+^ K 
can be formed by combining this result with the 
result for the semileptonic amplitude in ref. [7j. 
Since high precision is the ultimate goal, we will 
not quote a result in this early stage of the error 
analysis. 
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Figure 4. The D s leptonic decay amplitude as a 
function of sea quark mass. 



clusters deployed at Fermilab under the Doe Sci- 
DAC Program. 

REFERENCES 

1. R. A. Bricrc et al, CLNS-01-1742 

2. D. Besson et al. [CLEO Collabora- 
tion], Phys. Rev. D 68, 032002 (2003) 
|arXiv:hep^e"x /0305100|. 

3. C. T. Davies et al. [RT QCD Collaboration], 
|arXiv:hep-lat /03040'04l 

4. C. W. Bernard et al, Phys. Rev . D 64, 054506 

(2001) [arXiv:hep-lat /0104002| . 

5. Aida El-Khadra and Matthew Nobes, private 
communication . 

6. J. Harada, S. Hashimoto, A. S. Kronfeld and 
T. Qnogi, Phys. Rev. D 67, 014503 (2003) 
[arXiv:hep-lat /0208004| . 

7. M. Okamoto et aL, |arXiv:hep-la t/0309107| 

8. S. M. Ryan, Nucl. Phys. Proc. Suppl. 106, 86 

(2002) |arXiv:hep-lat/0111010| . 

9. M. Wingate, C. Davies, A. Gray, 
E. Gulez, G. P. Lepage and J. Shigemitsu, 
arXiv:hep-lat/0309092. 



